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Abstract 



Fernandez Comment [1] on our pseudo-perturbative shifted-^ expansion tech- 
nique [2,3] is either unfounded or ambiguous. 
In his comment [1] on our pseudo-perturbation shifted - I expansion technique (PSLET) 
[2,3], Fernandez strived to prove that (I) PSLET is just a version of SLNT, (II) it is not true 
that PSLET enables one to obtain more perturbation corrections than SLNT, and (III) it 
seems that SLNT ( and, consequently, also PSLET) is divergent. 

We explain below why we believe criticisms (I) and (II) to be unfounded and criticism 
(III) to be ambiguous. 

• Our statement " the difficulty of calculating higher-order corrections in SLNT through 
Rayleigh-Schrodinger perturbation theory (RSPT) results in a loss of accuracy" is clear and 
need not be misleading. We refer to the comprehensive, historical, account in the work of 
Imbo et al [4], indicating the actual novelty of SLNT ( which could handle, via RSPT, only 
the first four terms of the energy series). Fernandez and co-workers ( in [6-8] of [1]) have 
used the hypervirial perturbation method (HPM) to calculate higher-order corrections in 
SLNT. Therefore one would call their method HPM-SLNT, or, at least. Modified - SLNT 
(as they themselves named it) and not SLNT. 

• We did not claim that PSLET is completely different from SLNT [4] ( c.f. our comment 
following equation (31) in [3]). At the top of page 3063 in [5] we commented on the higher 
accuracy of the Fernandez HPM-SLNT method (although we had reservations about the 
order-dependent shift approach to the Klein-Gordon and Dirac equations). 

• Fernandez derived relations between a and /3, k and I, ■ ■ -etc. However, that work 
simply illustrates part of the message which we tried to deliver to readers, i.e., SLNT is not 
an expansion in large-N but, in effect, an expansion in large-£ ( c.f. C M Bender et al [6]); 
hence we preferred the abbreviation PSLET . 

• It is true, of course, that our conclusions in [3] about numerical accuracy referred 
to calculations for state wavefunctions with at most one node. However, the comment by 
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Fernandez that we are unable to apply our method to wavefunctions with more than one node 
is unjustified, since we have in fact given results for such functions in the tables of [2,7,8]. 
Below we also report PSLET results for the truncated Coulomb potential V{r) = —l/{r + a) 
with a — 10 for wavefunctions with several nodes. We show the sum of the first twenty terms 
of the energy series, £'20, and list the corresponding Fade approximants. The orders at which 
the energy series and Fade approximants stabilize are also shown. 

stabihty starts stabihty starts 

State £'20 Fade 

from from 

4s - 0.011638 E12 - 0.011638 £[3,3] 

6s - 0.006795 £12 - 0.0067958 £[7,8] 

7s - 0.005443 £12 - 0.0054438 £[7, 8] 

9s - 0.003721 £13 - 0.003722 £[7, 8] 

lis - 0.002705 £14 - 0.0027068 £[8,8] 

• Fernandez is unjustified in asserting that FSLET is based on logarithmic perturbation 
theory (LFT) (c.f. Appendix A in [4] and the references cited therein on LFT). FSLET is 
simply an algebraic recursion method which leads to exactly solvable recursion relations ( 
based on the uniqueness of power series representations, c.f. [9]). 

• It is not universally true that HFM-SLNT and consequently FSLET are diver- 
gent. Both techniques are based on asymptotic series expansions and one would expect 
to get asymptotically divergent or asymptotically convergent results (c.f. our analysis in 
ref.[7,10,ll]). To illustrate this statement with some persuasive evidence we consider the 
truncated Coulomb potential with a = 10, for wave functions with 10 nodes at £ = 1, 3, 5, 15. 
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Obviously, the trends of convergence are very well marked. In general the energy series 
of SLNT, HPM-SLNT, and PSLET are oscillatory ( a signal of, at least, asymptotic con- 
vergence) and one would, as a remedy, use an order-dependent shift ( as in HPM-SLNT) or 
Pade approximants ( as in PSLET) to obtain results with satisfactory accuracy. 

We agree with Fernandez about the unfavorable case { a — 0.1, £ — u — 0). Here the 
energy series appears to be asymptotically divergent. However, this should be attributed 
mainly to the nature of the truncated Coulomb potential and to the irrational value of 
a. One should notice that this particular potential gives contributions to the higher-order 
corrections of the energy series through its non- vanishing higher-order derivatives. This will 
lead to accumulated rounding-off errors which, in turn, can yield unreliable results from the 
higher-order corrections. 

We beheve that the points made above have satisfactorly answered the criticisms (I) to 
(III) of Fernandez. 
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